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Abstract. In this paper it is shown that the location of the characteristic overshoot 
of the Richardson's annular effect changes with the Kinematic Reynolds number 
CO* = urQ /v in the range of frequencies within the laminar regime. Prom the study of 
the Richardson's overshoot at different times it was identified the existence of aparent 
transverse damped waves similar to those ones observed in the famous Stokes second 
problem, the physical analysis of this waves was used for the establishment of a semi- 
empirical law that gives the functional relation of the mean overshoot maximum with 
the kinematic reynolds number, say Bo{oJ*) = 2.28 + 0.51w*~^^^. Finally a transient 
solution was constructed and verified asymptotically for large times, and the tipical 
time for which the transient solution resembles the steady oscillating one was identified 
to be dependent of the viscosity of the fluid and of the radius of the pipe. 
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1. Introduction 

Nowadays there is a little bunch of analytical solutions on transient flows, tipically this 
kind of solutions doesn't satisfy the initial condition of rest flow at t=0, including those 
ones related to pipe flows, specifically consider the well known velocity field of circular 
pipe fiow due to an oscillating pressure gradient. 

Is clear that circular pipe fiow motivates the use of cylindrical coordinates for the 
Navier Stokes equations. Now two fundamental hypotheses will be taken, first a regime 
in which the angular component of velocity vanishes = 0, and second a fiow that 
doesn't change in the axial component, say z; this means that the remaining components 
of velocity are given with the following dependences: 

Ur^Ur(r,t) 

Uz = Uz{r, t) 

The continuity equation in cylindrical coordinates under the two previous 
hypotheses gives: 

~{prur)^0 (3) 

Where p is the fiuid constant density. Now the no-slip condition for the radial 
component together with (3) show that the radial velocity must vanish everywhere 
inside the pipe, giving a one-directional fiow given by the axial component Uz{r, t). Now 
the Navier Stokes equations in cylindrical coordinates reduces to the simplified axial 
equation: 

fdu;, du^ duA dP d f duA d^uA , 

All the derivatives in z give zero using the second assumption of axial independence, 
and the remaining non-linear term in (4) vanish due to the fact that the continuity 
equation (3) gives = 0; finally writting u^ '■— u the axial momentum equation is: 

du dP ( d'^u ldu\ . . 

Pm^-Tz^^[d7-^-rd-r)- 

The most natural way to solve the steady oscillating state of (5), for which the flow 
is already oscillating with the same frequency of the external pressure gradient, is by 
introducing an harmonic velocity fleld: 

u{r,t) = R{r)e''^K (6) 
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By using the no-slip condition at pipe wall u{rQ, t) = 0, where vq is the pipe radius, is 
easy to check that the variable separation assumed above leads to the following velocity 
field [I]: 

K ■ , ( MrJ-iu/u) \ 
uir, t) = -e-* 1 - . (7) 

However this fiow doesn't satisfy the initial condition of rest state at t = just 
as was mentioned before; evidently only in the case when the pressure gradient starts 
with a value equals to zero the initial condition is satisfied, this is equivalent to have a 
sinusoidal gradient like sin{ut), but it is desirable to have a solution that starts from 
rest for a general harmonic pressure gradient as in (1). 

2. Richardson's Annular effect 

The particular oscillating fiow described by (7) exhibits a quite interesting behavior 
known as the Richardson annular effect yy, which consists of the presence of a velocity 
overshoot near to the wall, this can be observed for high frequencies of the pressure 
gradient in pipe fiow experiments; Sexl [3] observed theoretically this overshoot by 
expanding the solution (7) for large frecuencies, thereby he achieved the following 
asymptotic fiuid fiow: 



u{r*,t) 4 . , e 



-B 



sin{ut) -^sin{ut — B) 



(8) 



Where u* = utq/u is the kinematic Reynolds number of the problem, r* = t/tq is 
a dimensionless radial coordinate, Umax = Kro/ii' is the centerhne velocity for steady 
Poiseuille fiow with pressure gradient -pK, and B is defined in terms of ro and u* as 
follows: 




B = {l-r*)J-. (9) 



A plot of the mean square velocity derived from the asymptotic fiow (8) is given in 
figure 1 from which is easy to see the overshoot for B near to zero, that means that the 
overshoot occurs near to the wall of the pipe. 

The mean squared velocity was used instead of the mean velocity because due to 
the oscillating character of the last one the mean velocity gives u = 0. 

3. Transverse Damped waves in oscillating flow 



A plot of six velocity profiles given by (8) and equally time spaced over a period of 
oscillation for a kinematic Reynolds number co* = 500 is given in figure 2. 
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B= (1-r') V'w' / 2 

Figure 1. Plot of the Near Wall velocity overshoot for the specific case uj* — 500 



0.002 
0.001 

U 0,000 

-0.001 

-0.002 











t=T/6 






t=T/3 












t=2T/3 






t=ST/6 































O.S 0.6 0.7 0.8 0.9 1.0 



Figure 2. Equally time spaced velocity profiles over a period of oscillation 



It clearlly shows that the location of the instantaneous maximum of the overshoot 
changes with time. To observe with detail the location of maximum, the derivative of 
the velocity field in (8) was calculated: 

u-ir'A ^ A { _ B)\ . (10) 



Here the vanishing point of the derivative given by u'{r'^^^,t) = 0, indicates the 
location of the maximum. Plots of this function for the same times used in figure 2 for 
the velocity profiles are shown in figure 3 [|j They suggest that the overshoot maximum 
moves somehow periodically with time, being the period of the motion the half-period 
of the pressure gradient tt/u. 

During a half period, as time increases the value of r^^^ decreases, that means the 
the maximum of the overshoot is moving from the wall to the centerline of the pipe; 
then after the half-period the maximum of the overshoot appears again near to the wall 
of the pipe. 

I The crossings with the horizontal axis correspond to the solutions of u'{r*,t) — 
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It seems that it returns to its original position after a half-period of the pressure 
gradient (see intersection of continous and dashed lines of the same color in figure 3). 
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Figure 3. Equally time spaced plots of u'{r*) over a semi-period of oscillation. 



Indeed the location of the overshoot maximum rj^^^. is quite dificult to define in 
this case since the equation that define it as a function of time is given by a trully 
complicated non linear relation that arises from the condition u'{r^^^,t) = which 
turns out to be: 

-1 



tan{u)t — B) 




This equation has multiple inverse functions because for a fixed time there exist 
infinite locations where u'{r*,t) vanishes; it can be seen that each location and therefore 
each root of (11) moves along the r* axis towards zero as times passes, eventually one 
of this roots may take a value smaller than Tq = 1 which is the limit for this roots to 
have physical significance. Then what is really happening is that as time move forward 
an initial overshoot moves radially from the boundary towards the centerline of the 
pipe; as long as this overshoot moves it decays and become negligible, at the same time 
another new overshoot appears in the wall of the pipe and follow the same history of 
the previous one and so on. This analysis suggest a motion that behaves like some kind 
of wave moving radially to the centerline of the pipe. 

Actually this behavior corresponds to transverse decaying waves in the pipe. The 
word "waves" in this case is quite dangerous because this a problem of difussion instead 
of a wave propagation one, so these are not real waves in the physical sense, they just 
show that the diffusion effects inside the pipe propagate in a oscillating wa}|H however 
in what follows, for the sake of simplicity, we will use the term wave to refer to this 
phenomenon. 

§ Here it was followed the interpretation of Kundu & Cohen [5], where they discuss the oscillating 
plate flow that resembles the same kind of wave appearance. 
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The existence of these waves was noted by Stokes for the so-called oscillating plate 
flow, where he proposed a depth of penetration for them as [6]: 

Hence the amplitude of these waves is damped by a exponential factor e~^^^, this 
can be seen for this case of pipe flow from the velocity field (8) where we see a exponential 
contribution e~*-^~^ )\/^*/'^ , which gives always positive values of (1 — r*) inside the pipe, 
this can be identified with the radial distance traveled by the waves from the wall and 
this leads for our case to the same length defined by Stokes. 

Now it can be analized the behavior of the mean location of the maximum overshoot 
by computing the derivative of the mean squared velocity u'^; so the mean location can 
be founded from: 

d ( 2e-^ , , e-2^\ , , 

1 - -^cos{B) + — - = 0. (13) 



dr* 

Let Bq{ijj*) be the solution of this equation as a function of the Kinetic Reynolds 
number; this equation is non-linear so an analytical solution for the unknown Bq{u:*) 
cannot be given. 

Now the depth of penetration defined in (12) changes like 0{ijJ^^^'^) and it is 
plausible to assume that the mean location of the overshoot given by (13) is directly 
related with the depth of penetration, so in a tentative way it can be proposed a 
functional form for Bq as follows: 

i?o(a;*) = A+^. (14) 



Numerical solutions of (13) for different frequencies were obtained (see dots in figure 
4), also a non-linear least squares analisis using the function (14) was applied to the 
numerical solutions of (13), giving the parameters A = 2.28 and B = 0.51, the fitted 
function (dashed line) and the computed data matches quite good as shown in figure 4, 
the correlation coeficient of the fitted curve gives = 0.99998. 



4. Unsteady decaying State 

To obtain a solution of (5) that satisfies the initial condition of fluid flow at rest 
when t=0, is necessary to assume a more general ansatz, motivated by the form of 
the homogeneous solution given when dP/dx = [4], thereby it is proposed a solution 
of the form: 

oo 

n(r*,t) = ^a„(t)Jo(j„r*) (15) 

n=l 

Where Jo{x) is a Bessel function of order zero of the first kind and j„ is the nth 
zero of Jo{x). Now introducing (15) in (5) gives: 

E (a'nit) + ^a„(t)') JoUnT*) = i^e^"*, (16) 

n=l \ fo J 
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Figure 4. Mean overshoot maximum as a function of the kinetic Reynolds number, 
the dots represent the numerical solutions of (13) for different values of w* and the 
dashed line is the function (14) adjusted by the non-linear regression. 



This turns out to be the Fourier-Bessel expansion of Ke'^'^^, hence the coeficients of 
the expansion are given by [5]: 



<(t) 



'0 



a„(t) = Ke 



lo xjQ[jnX)dx 



(17) 



Solving the integrals above by using some well known identities gives: 

<it) + -i^an(t) = — — . (18) 

'0 Jn^iyjuj 

The last relation is a first order ordinary differential equation fo the unknown 
coeficients of the expansion a„(t); now the initial condition M(r*,0) = implies that 
a„(0) = for all n in the sum, therefore the solution of this initial value problem is 
quite easy and gives: 



a„(t) = 2K- 



Finally the desired solution is given by: 

e 



u[r ,t) = 2Kj2j-. — , -2/ 2v 7 / ■ ^ MJnr*). 



(19) 



(20) 



This solution gives rest flow at t = 0, and for t — )■ oo the real exponentials vanishes 
and we obtain a steady oscillating state: 



u 



{r\t) ~ 2Ke''''Yl 



(21) 



^1 {iUJ + Vjl/rl)jnJl{jn)' 

Rigorously, in order to have the steady oscillating solution (7) we should have the 
equality: 



E 



' ^ _ Jo{r^-iuj/u) 



(22) 
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Here we give numerical proof that indeed both expresions agree for all values of r 
where it were defined the following functions: 

n=l(^C^+#)jnJl(jn)' 



N I 




Jo(royQavV)^ 

Therefore the equality (22) is equivalent to the following limits: 

hm ^{D^{r*))=0 





lim 5>(A^fr*~ 



(23) 
(24) 

(25) 
(26) 



Where this limit must hold Vr* G [0,1]. Numerical evidence that support the 
convergence of these limits is shown in figures 5 and 6 for a frequency uj = 100, where 
it can be seen that as N increases the real and imaginary parts of (r*) approach to 
zero. 

Now notice that (20) have a superposition of decaying exponentials each one with 
a diferent tipical time of decaying depending on the corresponding j„, so we have a 
decreasing sucesion of decaying times: 



(27) 



All decaying contributions vanish faster than the lowest one of them given by ti , so 
one could consider that a tipical time for which we recover the steady oscillating state 
solution must be Ti. 
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Figure 5. Convergence of the Real part of the Difference funcion D^(r*) as A'' 
increases for u! = 100. 



In Figure 7 can be seen the comparison between the transient and the steady 
oscillating solutions for a time t = Ti, showing the agreement of both solutions at the 
tipical time for a frequency of w = 100. 



The Richardson effect and a transient solution of pipe flow 



9 



2.x 10" 

o 

.-( 

H 

-2.x 10-=*- 



• H=15 

* H=20 



*• ■ ■ ^ • m" ^ ^ ■ *• 



0.0 



0.2 



0.4 0.6 



OS 



in 



Figure 6. Convergence of the Imaginary part of the Difference funcion Di^{r*) as A'' 
increases for uj — 100. 
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Figure 7. Comparison of the Steady oscillating and transient solutions at the tipical 
time Ti for a frequency lj — 100. 



5. Conclusions 

The precise numerical solution of the equation that defines the maximum of the 
overshoot from the mean squared velocity field showed that the location of the maximum 
in terms of Bq changes with the kinematic reynolds number u* in the range of frequencies 
within the laminar regime. 

From the study of the Richardson's overshoot at different times it was identified the 
existence of some kind of transverse damped waves in circular pipe fiow similar to those 
ones observed by Stokes in the problem of the fiow generated by a oscillating infinite 
plate, this observation was used for the establishment of a semi-empirical law that gives 
the functional relation of the mean overshoot maximum with the kinematic reynolds 
number, say Bq{u*). 

A transient solution was constructed and verified asymptotically for large times, 
and the tipical time for which the transient solution resembles the steady oscillating one 
was identified to be dependent of the viscosity of the fiuid and of the radius of the pipe. 
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